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ABSTRACT 


This study is an analysis of high resolution measurements taken during the 
PATCHEX experiment. PATCHEX was designed to study shear related mixing in the 
upper thermocline. As a participant in the study, the R/V POINT SUR’s contribution 
was to survey an 8 kilometer square area centered on the R/P FLIP which was moored 
at 34°N127 °W. The survey which was conducted from 4-16 November 1986, 
gathered high resolution CTD, microconductivity and current velocity data to allow a 
detailed analysis of the small scale dynamics and thermohaline structure of the site. 
This thesis examined a subset of the R/V POINT SUR measurements to describe the 
vertical and horizontal thermohaline structure, surface fluxes and heat content of the 
upper 120 meters at the PATCHEX site. 

A two-layer mixed layer system existed during the first half of the study period 
due to the effects of precipitation, a 0.1 N/m* mean wind stress and strong diurnal 
heating. A 0.04 to 0.06 ppt salt deficit and 0.25°K temperature excess were measured 
in the surface layer which ranged from 10 to 35 meters depth. The depth of the 
seasonal mixed layer varied between 42 and 62 meters in response to the internal wave 
field in the upper thermocline. 

Significantly different horizontal gradients were analyzed in the time series of 
temperature and salinity at 10 meters and 40 meters. This motivated an attempt to 
map those gradients on to a uniform grid by objective analysis. The horizontal 
gradients at 10 meters were determined to decrease most rapidly toward the southeast 
corner of the survey area. Weaker gradients at 40 meters increased nominally from west 
to east This further illustrated the dccoupling of the upper and lower layers from 6-12 
October. 

Convective mixing from a 100 W/m mean daily maximum latent heat loss 
combined with a 75 percent insolation reduction erased the surface layer on 8 and 9 
October and continued cooling the entire layer. A net surface heat loss of 35+5 
MJ/m? (0.7+0.1 MJ/m? through the surface) balanced a 0.6 to 0.8 MJ/m? overall 
mixed layer heat loss. This confirmed that the surface flux was the dominant heat 


transfer mechanism from 8-12 October. 
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I. INTRODUCTION 


The ocean distributes its heat, salt and momentum through various dynamic and 
thermodynamic processes. These processes range over scales of seven or more orders of 
magnitude from basin wide, geostrophic circulations to centimeter or millimeter scales 
where mechanical energy 1s viscously dissipated. Studies such as Garrett and Munk 
(1979) have developed and refined models to investigate the spectral universality of an 
internal wave base state in the open ocean. Because of the ubiquitous nature of 
internal waves, their viscous dissipation is believed to play a kev role in maintaining the 
fluxes of heat, salt and momentum in the ocean. Gregg (1984) observed patchy, 
intermittent mixing events within the well stratified region of the upper thermocline. 
These events persisted for at least 3 hours at the same depth as a near-inertial internal 
wave packet. He noted that the average value of €, the dissipation rate, 1s relatively 
low in the upper thermocline and suggested that sustained shear events associated with 
long period internal waves might provide the most significant source of mechanical 
mixing energy in that region. Near-inertial internal waves, which are intimately tied to 
meteorological forcing, would thus play an important role in the production of 
significant mixing in the upper thermocline. 

PATCHEX emploved four ships to study turbulent mixing patches associated 
with significant shear events in the upper thermocline and to relate observations of 
those patches to the internal wave field. The goals of the experiment as elucidated by 
Gregg et al (1984) were: 

e to describe the nature and spatial distribution of the observed density and 
velocity fields; 

e to discern the spectral nature of the internal wave field in both the frequency and 
wavenumber domains; 

e to observe and quantify the transfer of kinetic energy through the relevant scales 
and attempt to quantify related changes in potential energy; 

e to examine the variability of high wavenumber shear associated with major 
mixing patches; and 

¢ to refine theories governing the generation and evolution of mixing patches, Le. 
to relate average mixing rates to sustained shear events within the highly 
stratified upper thermocline. 


Acknowledging the difficulty in separating the contributions of shear and 
diffusion from turbulent mixing, the principal investigators decided to examine a 
horizontally homogeneous site with a simple salinity structure, free from the effects of 
double diffusion or salt fingering. The limited range of the San Diego based research 
platform, FLIP, a primary participant and a towed vessel, was also a major concern. 
FLIP’s range limitations clearly created the possiblity of significant diffusive 
complications associated with mesoscale features near the California Current. A 
location at a latitude much higher than San Diego was eliminated since the thermocline 
would be deeper and less accessible as one moved northward. After studving satellite 
imagerv and a review of the ODEX and CALCOFI hydrographic surveys indicated 
that there were no logistically feasible sites with significantly less horizontal structure, a 
location near the MILDEX (Mixed Laver Dynamics Experiment) site, which was 
familiar to the principal investigators, was chosen. 

Since preliminary information could not ensure an ideal, uniform horizontal field 
with a stable salinity structure and minimal external forcing, the ex post facto 
assessment of site suitability is an important preliminary consideration to any 
turbulence study that might be undertaken. An examinination of the horizontal 
thermohaline structure of the site and a study of the balances of heat, salt and 
momentum in the mixed layer and the upper thermocline is a logical approach to that 
assessment. The balances of salt and momentum are not considered in this thesis which 
examines the elements of the heat balance which characterize the upper ocean forcing. 
It will also determine whether significant horizontal thermohaline structure exists in the 
upper 120 meters at the site. This study relies primarily upon the isolation and 
measurement of any significant mesoscale or sub-mesoscale temperature and salinity 
features. Equally important is the analysis of salinity and temperature profiles to 
locate and measure intrusive signatures in the upper thermocline. 

The overall heat balance equation of Large et al (1986) outlines four processes 


which combine to close the heat balance to an arbitrary depth, D. This system is: 


61(D) = Q, + R,(D)- A,(D) - A,(D). (1.1) 


10 


Werte Ss. of the equation represents the total heat content or 


oo , 
H(D) = pc, JT dz (1.2) 


and the R.H.S. includes: surface heating 


Q. = pe, Fat at (1.3) 
vertical diffusion 
NO) pew) pb (1.4) 
vertical advection 
0 
A,(D) = pc, Men 0,T dz OES) 
and horizontal advection 
Ay(D) = pc, (° ye vlondz. (1.6) 
-D 


For all of the formulations in the above system f is the density of seawater, c . is the 
specific heat of seawater, and T is the ambient water temperature at depth z at time t, 
to t. Also, the variables, V and w, represent the horizontal and vertical velocity 
respectively. 

It would be desirable to attempt a complete solution of the heat balance 
equation and perhaps an analagous salt flux system. The primary purpose of this 
thesis, however, is to analyze the surface forcing, heat content and horizontal structure 
to determine the extent to which mixing in the PATCHEX domain may be considered 
to be solely due to vertical shearing processes. These three contributors directly 


influence the transfer of momentum, salt and heat into the upper thermocline. The 


I] 


solution of equations 1.4, 1.5 and 1.6 is deferred since a complete closure of the heat 
budget is not attempted in this thesis. Absolute velocity data available from the R/P 
FLIP and temperature microstructure data from the R/V POINT SUR walle te 
employed in a subsequent study to estimate the remaining terms in the heat balance. 

Studies such as D’Asaro (1985) show that large scale surface forcing 1s 
intrinsically involved in setting up the near-inertial wave climate responsible for the 
hypothesized shear related mixing. He indicated that storm induced mixing extended 
5-10 meters below the mixed layer and was correlated with the phase of near-inertial 
currents. The total heat content and surface forcing terms are readily calculated and 
can be usefully compared irrespective of the other terms in the heat balance. These 
terms affect the temporal evolution of the temperature field and the depth of the 
seasonal mixed layer and hence, the temperature and momentum fluxes into the upper 
thermocline which directly affect mixing there. 

Horizontal temperature gradients can be estimated due to the spatial nature and 
high resolution of the R/V POINT SUR measurements. These data, which are 
described in Chapter II, provide a detailed analysis of the character and the spatial and 
temporal evolution of the horizontal thermohaline environment. This information 1s 
employed to directly determine if significant horizontal temperature gradients exist 
within the PATCHEX domain which may be related to mesoscale features. In 
summary, [he goals of this thesis are: 

® to examine the surface forcing and the heat content elements of the heat budget; 

¢ to map the horizontal temperature structure of the survey area and follow its 
evolution; 

€ to relate the evolution of the temperature field to the ambient forcing and the 
total heat content; and 

@ to observe the vertical coherence of significant horizontal structures through the 
mixed layer and into the upper thermocline. 

To effectively consider surface forcing, heat content and horizontal temperature 
Structure, the R/V POINT SUR data were analysed and displayed using a variety of 
methods and formats. Time series of surface forcing and cumulative heating were 
produced using a full stability corrected algorithm. The total heat content of the mixed 
layer was determined as a cumulative time series by solving equation 1.2. Horizontal 
temperature gradients and vertical coherence of those gradients were tracked by 


examining time series of temperature at constant depths of 10, 40 and 100 meters. 


we 


Where those gradients were significant, it was attempted to map them onto a grid of 
the entire PATCHEX domain with an objective analysis. All of these tools provide 
valuable information which sets the stage for the microstructure analvsis of the 
PATCHEX data. Only after these areas are developed in conjunction with the detailed 
oceanographic svnopsis which is provided in Chapter IV, can an accurate assessment 


of shear related mixing be undertaken. 
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HOP eESCRIPTIGN OFM EXPERIVEN? 


A. THE PARTICIPANTS 

Although this thesis concerns measurements made by NPS researchers aboard 
R/V POINT SUR, the concurrent measurements made by the participants in 
PATCHEX collectively illustrate the spatial characteristics of the site and are described 
below. The Scripps Institution of Oceanography research platform, FLIP, was the first 
participant to arrive on station at 34°N 128°W where she was moored for the 
duration of the experiment. The main instrument suite aboard FLIP consisted of four 
systems. The first was a pair of Seabird CTD’s for profiling to 600 meters every three 
minutes with 1 meter resolution for temperature and 3 meters for density. Four 
downward slanting incoherent doppler sonar transducers were mounted at 38 meters 
depth and oriented 55° down from the horizontal to provide 20 meter resolution 
velocity shear measurements from 75 to 1000 meters depth. To describe the large scale 
horizontal velocity field within the mixed layer, two horizontally mounted transducers 
similar to those employed in the slant array were used. A finer depiction of the near 
surface velocity field was cbtained from a coherent acoustic doppler velocity profiler 
which measured shear over a 50 meter range to | meter resolution. 

The next vessel to arrive was the R/V POINT SUR on 3 October. A detailed 
account of her role is summarized in the next section of this chapter. The R/V 
THOMPSON from the University of Washington arrived on 4 October to deploy three 
primary systems. A Richardson Number (RhiNo) neutrally buovant float equipped 
with velocity shear and stratification sensors was deployed approximately 4 kilometers 
from FLIP. A Benthic Acoustic Stress Sensor (BASS) on the float provided 0.03 
cm/sec resolution of the velocity field at 5 points over a 5 meter range. The float also 
carried an array of eight fasc response thermistors which were spaced vertically and 
loganthmically throughout the BASS field to provide fine scale temperature 
information over the same range. Temperature and velocity microstructure 
measurements were obtained by the next two systems: an advanced microstructure 
profiler (AMP) which provided repeated profiles of small scale shear and temperature 
gradient to a depth of 300 meters and a multi-scale profiler (MSP) which measured 
shear, temperature and velocity microstructure, and profiler orientation to 1000 meters 
depth. 


The last participant to arrive on scene was the USNS DE STEIGEUR with a 
team from the Naval Research Laboratory (NRL). The NRL instrumentation consisted 
piimamiaeol three smstems. Them unique contribution towtheseffortewas a towed 
thermistor chain which resolved 0.5 meter vertical sections of thermal structure. To 
accomplish this, the chain contained 180 thermistors spaced 0.5 meters apart and 8 
conductivity sensors spaced 2.5 meters apart. To provide finer resolution of features 
discovered by the thermistor chain, NRL also carried a High Resolution Vertical Array 
(HVR) mounted on a Sea Soar profiler. The HYR contained 30 thermistors spaced 0.1 
meters apart and 2 conductivity sensors. An RDI Doppler Current Profiler (DCP) on 
DE STEIGEUR provided velocity information. Its depth range was 200 meters, 


vertical resolution 8 meters and velocity resolution 2 cm/s averaged over | minute. 


B. MEASUREMENTS FROM k;V PT SUR 

The R/V POINT SUR departed Monterey, California shortly after noon on 2 
October 1986 and arrived to rendevous with R/P FLIP on the evening of 3 October. 
To provide oceanographic ground truth for the experiments to follow, a 20 kilometer 
CTD survey was made on 4 October. Eight casts to 1000 meters were made within the 
survey area using a Neil-Brown Instruments Systems (NBIS) CTD. Shortly thereafter, 
NPS researchers began conducting a detailed survey of an 8 kilometer square box 
centered on FLIP. This was accomplished using a procedure known as a toywo. In a 
towyo survey, a towed instrument package 1s alternately raised and lowered through 
the water column while the vessel is underway. This is done to acheive a three 
dimensional picture from a single sensor package. The towyo was deployed 
continuously throughout the experiment to gather detailed CTD, microconductivity 
and finescale shear profiles. Four systems were on the towyo package deployed from 
the POINT SUR. These systems included: / 

e a Coherent Acoustic Doppler Velocity Profiler (CADVP) capable of estimating 
3-component velocity profiles 5 times per second with 10 centimeter resolution 
over a 5 meter range. It employs a projector and four receivers as illustrated in 
eure 2.1; 

¢ a Cox style microconductivity sensor sampled at 1024 Hz; 

¢ an NBIS CTD sampled at 32 Hz; and 


¢ pitch, roll, heading and 3 axis acceleration sensors. 
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Figure 2.1 R/V Point Sur PATCHEX TOWYO Profiler. 


Due to high platform motion the value of the velocity data obtained by the coherent 
system requires further processing and evaluation to determine its resolution of the 
finescale shear field. A ship mounted Ametek Straza 4400 Acoustic Doppler Velocity 
Profiler (ADVP) concurrently obtained velocity profiles. This 300kHz, 4 beam 
“JANUS” geometry system sampled every 0.6 seconds to provide velocity information 
to 102 meters depth. Dead reckoning for shipboard navigation was also augmented by 
the ADVP which provided accurate estimates of the ship's velocity relative to a 
reference depth current. A Sperry gyro-compass, Magnavox MX-11S80 transit satnav 
receiver, Internav LC408 Loran-C receiver and Motorola Miniranger (providing 
ranging from FLIP) completed the navigation suite. Thirty second averages of standard 
meteorological and oceanographic parameters, which are summarized in the next 


section, were recorded throughout the cruise. 
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The primary sampling strategy covered a butterfly shaped pattern whose 
southeasterly and southwesterly diagonals intersected alongside FLIP and whose sides 
ran in a north-south direction. It was employed to provide a repeated, quasi-synoptic 3 
dimensional representation of the velocity and thermohaline character of the upper 
ocean over an 8 kilometer domain including radial intersections with LIP. These 
intersections allow decorrelation scales for both thermohaline and shear features to be 


determined by intercomparison of concurrent POINT SUR and FLIP measurements. 


Start/end 





Figure 2.2. PATCIHEX Survey Track, 05-13 October. 


This geometric pattern which ts depicted in Figure 2.2 was continuously surveyed from 
the beginning of the exercise on 4 October to 13 October. On 13 October, a triangle 
pattern which eliminated half of the original track was adopted to decrease the cycle 
time between patterns. Figure 2.3 depicts this pattern which necessarily decreased thie 
coverage area in order to improve the synopticity of the sampling process. Until 13 
October, most of the towyo profiles taken were between 5 and 120 meters although 
deep profiles to 300 meters and 15 meter constant depth runs were also periodically 


made. After 13 October, nearly all profiles were taken from 100 to 300 meters. During 
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Figure 2.3) PATCHEX Survey Track, 13-16 October. 


the latter period, measurements were also taken in a drifting mode for several hours at 
a time in conjunction with MSP drops from the R/V TITOMPSON. Late on I1 
October and into 12 October, 3 additional 1000 meter CTD casts were taken along a 
20 kilometer diagonal which bisected the PATCHEX domain. 

A large scale velocity survey was conducted on Ii October in which 12 
expendable current profilers (XCP’s) were dropped every 4 kilometers from the POINT 
SUR along an L-shaped track. These XCP’s provide arbitrarily referenced horizontal 
velocity profiles with [0 meter resolution or better from the surface to 1200 meters. 
lour of these probes were beyond shelf-life and failed. The results of that survey will 


be encorporated in future PATCIIEX analyses. 


C. DATA ACQUISITION AND PROCESSING 

The PATCIIEX data set was large an diverse. Over 5 gigabytes (billion bytes) of 
data were recorded on 108, 10 inch, 9 track tapes. To collect and monitor the data 
streain effectively, an integrated data logging system was developed. Prior ficld 
programs by NPS researchers aboard the POINT SUR used the SAIL (Serial ASCII 


Instrumentation Loop) Data Acquisition System, or SDAS, to monitor standard 
shipboard instrumentation. PATCHEX also employed SDAS but modified it to output 
its data as a synchronized common header to the towyo logging system. The 
parameters which were recorded by SDAS are listed in Table |. Towyo data have 


alreadv been described. 


TABLE 1 
R/V PT SUR STANDARD INSTRUMENTATION 


“A 


Data : Data 


Year ; Pumped Flow Rate 
Month and Day : Satnav Degrees Lat. 
Hour and Nlinute ; Sauiav Minutes Lat. 
Second Satnav Degrees Long. 
LORAN Degrees Lat. : Satnav Vlinutes Long. 
LORAN Minutes Lat. : Ship’s Heading 

LORAN Degrees Long. : Ship’s Speed 

LORAN Minutes Long. : Time 

Dewpoint Boom SST 

Air Temperature : Visible Insolation 
Relative Humidity : IR Insolation 

Pumped Sea Temperature : IR Pyranometer Temp. 1 
Pumped Conductivity IR Pyranometer Temp. 2 
Pumped Salinity ; Forward Velocity 

Open : Mliniranger-A Range 
Relative Wind Speed \liniranger-B Range 
Relative Wind Direction 


l. 
Z. 
3. 
4. 
> 
6. 
i. 
8. 
2. 





The data, whose acquisition is summarized in Figure 2.4 were recorded on 10 
inch magnetic tape for later processing with a Hewlett-Packard 9020 minicomputer. 
ADVP data, which had been provided to the towyo logging system by SDAS were fed 
to the program PAMSUR which produced velocity profiles relative to the ship. 
PROCCTD processed towyo data that had been acquired simultaneously with SDAS 
and ran in a seperate partition of the 9020. The program provided spectral corrections 
to conductivity data from the profiler’s NBIS CTD. Given the faster response of the 
conductivity sensor compared to the thermistor, which had a slower and velocity 
dependent response, it was possible to provide corrections which ultimately resulted in 
6 parts per million resolution in 0.25 meter bins. The final result of this process was 
that salinity signatures were Kept free of spiking to minimize false overturn signatures 
when density profiles were estimated. To calibrate the microconductivity sensor, CTD 


conductivity data was high pass filtered and compared with the low passed 
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microconductivity signal. The processed profile data with corresponding directory 
information were stored for each profile on a Bering 70 MB hard disk. Figure 2.5 
details the steps in the data processing procedure. 

Accurate reconstructions of the ship’s track, which was dead reckoned relative to 
FLIP with course, speed and radar positioning underway, was finahzed using 
NAVEDIT, a Basic routine run on the HP-9020. NAVEDIT, which was developed and 
applied by J. Stockel at NPS, employed all of the navigation information acquired by 
SDAS. It selected the most applicable and precise data among the variety of navigation 
information in the data set and tailored each profile’s position along a “best track”. 

The raw physical parameters on hard disk were accessed by OPPROF. OPPROF 
is a multi-function, post-processing program written in Basic for the HP-9020. It 1s 
capable of accessing specified series of profiles and can output graphics and time series 
of both raw and computed oceanographic parameters. Table 2 is a ment of these 


outputs. The program can also provide tabular output or a metafile graphics format for 


TABLE 2 
OPPROF PROCESSING OPTIONS 


Zz 
= 


Data io. Data 


Temperature U Shear 

Conductivity ‘ V Shear 

Salinity Speed Shear 

Sigma T ; Buoyancy Frequency 

Optical Transmission (n/a) : Gradient Richardson Number 
Flourescence (n,a) : at & S Plot 

Chi , Chi and Variance Plot 
Demeaned U : T.S Plot 

Demeaned V : Raw U 

Demeaned Speed 


ie 
Z. 
2: 
4, 
a: 
6. 
7 
8. 
9: 
10. 





plotting on the IBM 370 mainframe at NPS. The entire series of programs for the 
FIP-9020, with the exception of NAVEDIT and some utility routines in OPPROF, were 
developed and written by [T. Stanton at NPS. 

SDAS data from the acquisition process outhned in Figure 2.4 reside on both 10 
inch tape and 3.5 inch diskette media. The data on the diskettes were downloaded to 
mass storage volumes for access on the NPS mainframe. Surface forcing time series 
were then produced on the mainframe with a FORTRAN program, PATFLUX. The 
stability corrected flux algorithm in PATFLUX 1s a subroutine written by P. Guest at 
NPS. 
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Figure 2.5 PATCHEX Data Processing I’lowchart. 
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A. SYNOPTIC METEOROLOGICAL SITUATION 

Large er al (1986) observed that over 63% of the autumnal cooling of the sea 
surface during the 50 days of STREX (Storm Transfer and Response Experiment) was 
due to nine 1-2 day storm episodes and only 23% was from radiative forcing and 
entrainment. Davis et al (1981) had previously observed that the evolution and decay 
of the mixed layer is also primarily a result of storm activitv. Near gale conditions with 
significant precipitation were noted from 1-3 October, 1986 at the PATCHEX site. 
Conditions had begun weakening, however, with light intermittent precipitation and 
clearing when the R/V POINT SUR began surveying the area at 1000 local time on 4 
October. By 6 October, fair conditions developed and continued for the rest of the 
experiment. 

The high winds and precipitation observed at the outset of the experiment were 
not related to a typical transitory cyclone in the mid-latitudes. Rather, they resulted 
from instabilities associated with the gradient between an intense northeastwardlv 
moving 1034.5 mb subtropical high in the central eastern Pacific near 42N 142W and a 
strong inverted trough of low pressure whose 1004.8 mb center was located over 
southern Nevada. The trough extended northeastward through California and Oregon 
and out to sea as the low center moved slowly eastward and weakened. The evolution 
of the surface pressure and wind fields at the outset of the PATCHEX experiment is 
depicted in Figure 3.1 The top panel is the surface synopsis at 011200Z October 1986. 
It illustrates the 30 knot northerly winds that prevailed at 35N 127W when the R/V 
POINT SUR departed Monterey. These winds persisted at the PATCHEX site (marked 
by an X in the figure) for over 24 hours then gradually abated as the pressure gradient 
weakened. The synoptic meteorological situation at 041200Z shortly after the arrival of 
the R’V POINT SUR and 5 hours before she began surveying is depicted in the 
bottom panel of Figure 3.1. At that time the winds were still northerly but at onlv 10 
knots. The top panel of Figure 3.2 depicts the surface synoptic situation at 060000Z. 
This is roughly when this thesis begins, with the winds increasing to 16 to 18 knots by 
070000Z. These winds persisted for 24 hours then decreased to 10 to 12 Knots by 
O80000Z. This can be seen in the bottom half of Figure 3.2. 
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Figure 3.2 06 and 08 October synoptic situation at OOZ. 
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The wind remained near 12 knots through 9 October and is depicted as such in the 
101200Z map, the top panel in Figure 3.3. By 110000Z the wind began to weaken and 
despite some gustiness at 111200Z, the wind continued decreasing and had become 
light and variable by 121200Z. The bottom panel of Figure 3.3 illustrates the synoptic 
situation at this time. Note the strength and extent of the semi-permanent subtropical 
anticyclone which is the underlying cause of the prevailing conditions. It is the 
extreme atmospheric subsidence associated with the subtropical high that creates a 
large dewpoint depression and sets the stage for substantial latent heat release from the 
ocean surface. Detailed time series of the actual forcing of the marine boundary layer 
by the transfer of ambient mechanical and thermal energy are explored in the next 


section. 


B. SURFACE FLUX THEORY AND BULK PARAMETERIZATIONS 
The ocean is forced by many thermal and mechanical processes. In the absence 
of precipitation, local wind stress and heat fluxes are the mechanisms primarily 
responsible for short term variations in the diurnal mixed layer depth and the position 
of the seasonal thermocline. The rain which was observed prior to the arrival of the 
R/V POINT SUR had ceased shortly after she began surveying. To determine the 
atmospheric forcing, continuous measurements of bulk aerodynamic properties were 
taken onboard the R/V POINT SUR throughout the cruise. Since direct measurement 
of vertical gradients of velocity, temperature and humidity at scales small enough to 
infer the fluxes within the surface layer of the atmosphere are generally not attempted 
at sea, thev were estimated by employing Monin-Obukhov similarity theory, which for 
this development, uses the notation of Panofsky and Dutton (1984). 
The small scale gradients in the atmosphere are a function of wind stress, heat flux 
and surface roughness alone. The effects of thermal and mechanical forcing on 
turbulence can be combined in theory to produce a scale velocity, ux, which 1s the 


Square root of the wind stress (t) over density (p), ie. 


and a scale length 


L = L(us, Qo) (3.21) 


oe 


which is called the “Obukhov” length. In turn, these fundamental quantities may be 
applied to develop temperature and humidity scales which allow the surface fluxes to 
be estimated: 


Altering the form of our velocity scale expression as stated in equation 3.1 to 


Us? = tip = <-uw> (3.3) 


fundamentally defines the wind stress. Monin-Obukhov scaling also develops a similar 


temperature flux formulation 
ligt > = < -Q'w’ > ; (3.4) 


where @’ the potential temperature, is the temperature that a parcel of air would have if 
it were displaced adiabatically from its present level to the 1000 mb level. There is also 


an expression for the specific humidity flux 


Usdx = <-q'w'>. (3.9 


A rigorous derivation of the Obukhov length scale is obtained by manipulating 
the turbulent kinetic energy equation. The solution of a scaling problem, however, is 
not achevieved explicitly. In practice, the velocity scale is determined first by initially 
expressing the wind shear as a product of scaling parameters and a “dimensionless 
Shear” such that 


du/0Z = @ Ux/kz (3.6) 
where k = 0.4 is von Karman’s constant and @, represents the dimensionless wind 
Shear. An accepted convention as expounded by Panofsky and Dutton (1984) states 
that 

yp, = 1-166, 25§ <0 (3.7) 


under unstable surface conditions, and 


p_, = 1- 7, E> 0 (3.8) 
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for the stable case. In these expressions § = z/L. The constants 16 and 7 continue to 
be the subject considerable debate and research. Large and Pond (1982), which this 
study follows, suggests that for the stable case the constant is 7, and not 5, which ts a 
convention used by Panofsky and Dutton (1984). The constant in the unstable case 1s 
generally agreed to be 16. 

Integrating the expression 3.6 from the surface roughness height (Z)) to the 


measurement height (z) produces the diabatic profile relation 


ue = ku { In(z/Z9) - v,,(8) } (3.9) 


in which y_(§) can be approximated for the unstable case with the formula of Paulson 
(1970) 
w= In ((1+X?)/2 + ((1+X)/2)7} - 2 atan X + 2/2, (3.10) 


m 


where the variable, X, is defined as 
X = Ay (l= 162) (Sabily) 


For the stable case, a much simpler relation exists, i.e. 


wh = -76. (3.12) 


Three methods are commonly used to evaluate scaling arguements. These 
include the eddy correlation method, the dissipation method and the bulk aerodynamic 
method which is employed in this study. Since @, = @_(§), as noted in equation 3.7 
and € cannot be found explicitly, an iterative technique is used to refine its value while 
solving equation 3.9 for ux. The process begins with an initial determination of the 


roughness scale (Zp) using a neutral drag coefficient, CDN}, ie. 


Zy = z exp(-k//CDN j 9). (3.13) 
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The drag coefficient of Large and Pond (1982) is used. It is a wind dependent quantity 


whose form in lighter winds 1s 
CDNy9 = 1.14x 10°, U < 10m/s. (3.14) 
In heavier winds equation 3.14 becomes 
CDNjo = (0.49 + 0.065U) x 10°, 10 < U < 25 mis. (3.15) 


Applying the method of Large and Pond (1982) PATFLUX performs an iterative 
procedure on Tx and Qs to provide a stability correction for §. Iteration stops when 
the algorithm fails to improve us by more than 5%. The wind stress, or momentum 
flux was obtained by applying equations 3.9 and 3.3. The net radiative flux is defined 


as 


Qy = (1-€)Qg + Qy - Qp - (Qe + Qzy) (3.16) 


where @ represents the albedo or average reflectivity of the ocean surface estimated at 
15% for the PATCHEX cruise. This value is based on a graphic global summary of 
albedo contained in Gill (1982). Qe and QL, incident solar visible and longwave 
radiation Were measured by two shipboard pyranometers. These instruments were 
mounted on gimballs atop the bridge cabin to ensure that the sensors continuously 
looked overhead. Infa-red radiation emitted from the sea was estimated using 


Stephan’s radiation law for black bodies 

Qp = eoT.? (3.17) 
where the emissivity, € = 1.0, to properly account for reflected longwave (Large et.al., 
1986). Stephan’s constant or © = 5.67 xX 10.9 W m* K~4. T, is the sea surface 


temperature or SST. The evaporative or latent heat flux was evaluated from the 


expression 


QE = PL usQu. GAS 
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Sensible heating was determined using 


QT] = “Pc us Ts. ovr} 


In equations 3.18 and 3.19, p is the density of seawater. The temperature is expressed 


in absolute degrees, i.e. 
ele =a Gy 27 aelo (3.20) 


C. ESTIMATING THE FLUXES 

Accurate measurement of the meteorological and radiative forcing at sea is 
particularly challenging. The heat production and storage, stabilitv, and motion of the 
platform all contribute to this difficulty. Monin-Obukhov scaling, the standard 
approach to providing flux estimates, is very sensitive to errors in the four primary 
input parameters: air temperature, dew point, sea surface temperature and wind speed. 
The more extreme criticisms of current application of scaling considerations, such as 
Blanc (1987), belabor worse case senarios where measurement errors produce flux 
errors approaching 2 to 3 orders of magnitude. However, until there is quantum change 
in our understanding of air-sea interactions and in the technology employed at the air- 
sea interface, the present approach provides one of the soundest estimates of the fluxes 
of heat and moisture (Davidson, personal communication). 

The flux of latent heat from the sea surface is best determined by measuring the 
surface temperature in situ or by direct pyranometric measurement. Both methods have 
their advantages and drawbacks. It is difficult to measure the skin temperature in situ 
Without penetrating the interface. On the other hand, pyronometeric readings also 
contain longwave energy from the sun and clouds which is reflected off the sea surface. 
Pvranometric SST measurements were not employed in PATCHEX. The R/V POINT 
SUR’s stern boom trailed a floating platinum resistance thermometer having a 0.001°C 
overall temperature resolution. It was expected to provide the most accurate SST 
measurements. Pumped seawater from several meters below the water line was also 
measured in an onboard sea chest to effectively provide a “bucket temperature” as a 
back-up to the boom temperature. Both of the measurements described are superior to 
engine injection temperatures used in some older studies, which fail to reflect the true 
SST due to heat contamination from the engine room and due to the subsurface depth 
of the intake. 


oI 


It was essential to employ the boom SST data because of the surface 
Stratification which existed after the period of precipitation and during intense daytime 
solar heating. The preliminary audit of the data set had determined that the difference 
between the boom temperarure and the CTD temperature was relatively constant and 
reflected an offset in the boom measurements. Both types of POINT SUR’s SST 
measurements Were compared against toywo temperature measurements in the upper 5 


meters of the ocean during convective periods when the entire mixed layer was 


VNB EE. 
SS PMIEASUREMENT DIEPEREN GES, 3 OGTOCEK 3-5 Acy1 


06 Oct 08 Oct 10 Oct bP 2Ocs 


Seachest 0.077 0.075 0.040 0.032 
Boom 0.208 0.161 0.157 0.162 


Mean AT: Seachest 0.055 Boom 0.172 





isothermal. Table 3 is a summary of the differences between the SST and CTD 
measurements from approximately 3-5 AM local time on 6, 8, 10 and 12 October. The 
boom temperature difference is approximately 3 times the magnitude of the sea chest 
difference. Figure 3.4 details the range of temperature differences that were observed 
for the two sensors. Those ranges are 0.024 and 0.026 for the boom and seachest 
respectively. The result of the towyo comparison determined an offset to be applied to 
the time series of the boom SST of -0.172 °C. 

Air temperature and wind speed data were not corrected, though the dewpoint 
was estimated from psychrometer readings for a 24 hour period while the dewpoint 
sensor was down. The raw forcing data were then used to produce time series of the 
forcing, the fluxes and the cumulative fluxes for the study. Hourly averages of the 
meteorological forcing are summarized for the entire cruise in Figure 3.5. This thesis 
covers Julian days 279 through 285 which 1s the period from 6-12 October 1986. 

A general description of the winds was provided in the first section of this 
chapter. The detailed time series shows that although there was some variability, the 
wind was generally fresh and remained between 6-8 m/s for most of this study. This 


was a critical factor in the maintenance of daytime convection throughout the oceanic 
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Figure 3.4 SST Measurement Differences. 


fossil muxed layer. Toward the end of the experiment, as subsidence associated with the 
prevailing high pressure peaks, the time series evidences a sharp decrease in the wind. 
‘\ steadily increasing dewpoint depression can be seen in the time serics of hourly mean 
air temperature and dewpoint. By the end of the study period, it had reached 8°C. 
Liaison with FLIP (Pinkcl, personal communication) confirmed that less than 50% 
relative humidity persisted over the last half of our measurement period. A peak in 
SST occurred on October 6 and is depicted in the timescrics of boom and sea chest sea 
surface temperature. The peak corresponds to the period of maximum surface rainfall 
accumulation. Note that on previous days, night-time convection had mixcd the 
surface water with the seasonal mixed layer. By the end of 6 October, this mixing had 
erased all evidence of a surface layer. Smaller maxima in SST and a gradual cooling 
trend were then observed in tlic remainder of the timeseries. 

The combination of low relative humidity and moderate winds removed a large 
amount of heat from the sca surface. Although the mean SST depicted in Figure 3.5 
only reficcts a 0.2-0.3 °C change from the evening of 6 October through the end of 12 
October, this represents a significant heat release. Figure 3.6 is a timesceries of the heat 
and momentum fluxes for the duration of the PATCIIEX cruise. Since our period of 


interest is from Julian day 279 to 286, note the timescries of latent and sensible heat 
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Figure 3.5 PATCHEX Meteorological Forcing. 


loss depicted in the two bottom panels of Figure 3.6. Two peaks of 125 W/m? of latent 
heat release can be seen. These maxima occurred tn conjunction with increased wind 
stress episodes on 6 and 7 October. The mean latent heat loss until [1 October was 
approximately 100 W/m*. Corresponding sensible heat losses on the order of 75 W/m? 
can be seen tn the timeseries. The winds began to decrease on Julian day 285, or 12 
October and by day 286 sensible and latent heat loss had dropped to 50 W/m? and 20 
W/m? respectively. 
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Figure 3.6 PATCIIEX Surface Fluxes. 


Combined solar longwave and shortwave cnergy, minus the output of longwave 
energy from the sea surface, is depicted in the third panel of Tigure 3.6. The typical 
diurnal cycle shown ranged from approximately -200 W/m? overnight to 500-700 W/in? 
during peak daytime heating periods. This cycle was interrupted on 9 and 10 October 
when cloudy conditions limited insolation. In the next chapter, we shall observe that 
this period of retarded surface heating was instrumental in sustaining and intensifying 


the cooling within the scasonal mixed layer which erased all traces of a surface layer 
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(even during the day). The time series of net flux, Q,, depicts the overall heating and 
cooling which was estimated at the air-sea interface. It is clear that solar radiation 1s 
the largest contributor during the daytime. It should be noted that night-time losses of 


net heat were increased 3 or 4 times due to the effects of sensible and latent heating 


after 8 October (Julian day 281). 
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Figure 3.7 PATCHEX Integrated Surface Iluxes. 
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The surface fluxes were integrated for the study period for later comparison with 
the total heat content of the fossil mixed layer. Along with the hourly mean wind 
stress, these quantities are shown in Figure 3.7. Following a & 5 MJ;m? diurnal cycle 
through 8 October a -35 MJ/m* change in the net flux is observed through 11 October. 
This is followed by a net increase of approximately 5 MJ/m? on 12 October. The wind 
Stress and sensible fluxes associated with the subtropical anticyclone were clearly 
insufficient compensation for decreased insolation and continued latent heat release. 
The cumulative time series for Julian day 286 through 287 depict additional heating but 


are not included in this study. 


D. ERROR ANALYSIS OF THE SURFACE FLUXES 
The accuracies of the various instruments used to measure raw flux paramcters 


are listed in Table 4. They are based upon pre-cruise calibrations for the pyranometers 


TABLE 4 
ACCURACY OF FLUX SENSORS 


Sensor Accuracy 


Air Temperature 0.2 nC 
Dew Point Temperature O72 

Sea Surface Temperature 0.01 C 
Wind Speed lms 
IR Pyranometer eo 
Visible Pyroheliometer 0.5% 





and both pre-cruise and post cruise calibrations for the temperature sensors. The 
windspeed error is based upon the anemometer specifications from the manufacturer. 
The values in Table 4 were employed along with environmental considerations and an 
audit of the time series of raw parameters to estimate the errors listed in Table 5. 

The 30 second average values, recorded by SDAS, provided some smoothing of 
the raw data from the sensors. For example, although the time series of wind stress 
appeared noisier during the first days of the experiment, the variability was within 
range of the ambient environmental conditions and not attributable to instrument 
error. Low sea state also minimized pitch and roll pumping effects on the anemometer. 
The error in t is based solely on instrument error as a percentage of the total wind 


stress. 


a1 


Along with the wind stress, Table 5 also contains | percent error for the fluxes of 
incoming and outgoing lor.gwave energy. Incoming solar infared energy error was 
based upon the | percent pre-cruise calibration of the pyranometer. Surface emissivity 
or back infared radiation, was determined strictly as a function of the sea surface 
temperature where the boom was recalibrated to within 0.1°C following the cruise. 
Incident shortwave solar energy error was estimated at 2 percent, combining the 0.5 
percent advertised error of the pyroheliometer with an estimate of the effect of mast 
shadowing. This effect was estimated as a percentage of outliers in the time series of 


visible solar flux during clear sky conditions. 


TABLE 5 
ACCURACY OF SURFACE FLUX ESTIMATES 


Percent Error 





Sensible heat flux was estimated based upon the combined effects of the error in 
air temperature, SST and wind stress as a percentage of the total effect. This accounted 
for 3 percent error. The largest error term, Qr, the latent heat, had the addition of the 
dewpoint error which raised its error to 5 percent. The total error of all the fluxes when 
combined is 14 percent. This equates to a cumulative flux error of +5 MJ/m? for the 


35 MJ/m? which was lost at the sea surface during this study. 
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IV. OCEANOGRAPHIC SYNOPSIS 


A. VERTICAL CHARACTER OF THE PATCHEX SITE 

Chapter II described the spatial and temporal nature of the R/V POINT SUR 
measurements which allow the thermohaline character of the mixed laver and the upper 
thermocline of PATCHEX site to be determined. This high resolution dataset provides 
a detailed oceanographic context for the pursuit of the goals of this study and for the 
analysis of the microstructure and internal wave measurements made by all of the 
participants in PATCHEX. [In this study, the decision to explore the different areas 
described in Chapter I was based directly on an analysis of the thermohaline 
environment. The first phase of that analysis concentrated on the evolution of salinity 
and temperature profiles within the seasonal mixed layer from 6-12 October 1986. 

The preliminary analysis first shows that from 6-8 October the mixed layer was 
undergoing incomplete daytime mixing due to precipitative stratification and weak 
surface fluxes. Convective mixing of the entire layer was observed through the night 
by the absence of jumps in temperature and salinity to the resolution of the CTD 
(0.001°C and 0.2 millisiemens respectively). Figure 4.1 illustrates typical daytime 
salinity and temperature profiles observed during the first two days of this study. 
Temperature steps where AT < 0.05°C was observed generally coincided with an 
active mixing layer. These data were not, however, as sensitive a test for well mixed 
conditions as was the salinity signature. To determine whether or not thorough mixing 
existed, salinity was used as the primary tracer. A survey of all of the profiles in the 
experiment indicate that if AS < 0.04 parts per thousand, the entire seasonal mixed 
layer was overturning. 

The salinity signature clearly indicates that active mixing is limited to the region 
above 10 meters in Figure 4.1. Typically these salinity ledges ranged from 10 to 35 
meters in depth during the day. To ensure that salinity features such as that seen in 
Figure 4.1 are not dismised as sampling artifacts, it is reiterated that 0.5 meter vertical 
salinity resolution of 0.003 parts per thousand was obtained in the pre-processing 
phase. It should be noted that these salinity jumps were frequently horizontally 
coherent features which persisted from profile to profile throughout the domain over 


time scales spanning several hours. Salinity ledges such as the AS of 0.06 ppt near 10 
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Figure 4.1 PATCIHIEX Daytime Salinity and Temperature, 6 October. 


meters in Figure 4.1 were present during periods of strong surface heating. Otherwise, 
the step was generally observed near 25-35 meters combined with the AS normally seen 
at that depth (0.03ppt at 35 meters in the example) forming a single, stronger step at 
the deeper depth. It was common to observe completely active surface layers with a 
single AS as high as 0.16 ppt near 25-35 meters. These events occurred sporadically 
throughout the day in conjunction with the tenuous balance that existed with the 
overall forcing. An active daytime surface layer such as this was not observed after 8 
October when mixing from wind stress and strong latent heat loss kept entire layer well 
mixed. As seen in Chapter III, the wind stress quickly peaked on 8 October with 


corresponding sensible and latent heat losses to begin re-mixing the stratified surface 


40 


layer tnto the fossil mixed layer below. After 8 October, the mixed layer continuously 
exhibited an tsothermal and tsohaline character. 

The average depth of the base of the mixed layer ranged from 48 to 62 meters. 
The effects of tnternal waves were clearly a dominant mechanism in the vertical 
displacement of the base. Figure 4.2 depicts a typical Brunt-Vaisala profile from the 
experiment. The large stability signature characteristic of the uppermost boundary of 


the thermocline tllustrates the favorable climate for tnternal wave activity there due to 
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Figure 4.2 PATCHEX Buoyancy Profile, 6 October. 


the very strong stratification. Figure 4.3 is a time series of the depth of temperature 
isotherms for 12, 13, 14 and 16°C, based on averages of every twenty profiles from 


6-12 October. The average elapsed time between each profile was between 6 to 9 


4 


minutes. The crude low-pass filter produced by the averaging of profiles removes the 
noisy higher frequencies with periods less than 3 hours and shows the underlying near- 


inertial or semi-diurnal periodicity by inspection. Consecutive semi-diurnal and tnerttal 
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Figure 4.3 12, 13, (4 and 16°C Isotherms, 6-12 October. 


periods are denoted in Figure 4.3 by the symbols d, and f, respectively. The inertial 
period at the PATCIIEX latitude ts 20.84 hours. Dashed lines represent consecutive 
20.84 hour tncrements while alternating dashes and dots represent consecutive 12 hour 
increments. Inspection of the time series appears to indicate that neither an inertial nor 
semi-diurnal period dominates the upper thermocline from approximately 50 to 100 
meters during this study. For the purpose of this thesis, it is important to establish 
that below the mixed layer, temporal variations may overshadow the spatial variations 
that are investigated tn the next section. The effect of internal waves ts also reflected in 
the deeper, slightly stratified areas of the mixed layer which appears to contatn less 
horizontal coherence. 

To properly characterize conditions within the mixed layer, the complex 
thermohaline structure in the active upper thermocline cannot be tgnored. The activity 
near the top of the thermocline during PATCHIEX ts thlustrated in Figures 4.4 through 
4.7 at the end of this section. Temperature inversions on the order of 0.5°C, as seen in 


Figure 4.4, were common directly below the base of the mixed layer. These inversions 
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occurred throughout the study and were often as strong as 1°C. Figure 4.5 depicts the 
Brunt-Vaisala frequency for the same profile emphasizing the strong stratification 
which exists in the upper thermocline. Figure 4.6, a velocity profile based on the 
relative velocity information obtained by the Amatek Straza doppler system depicts a 
15 cm/sec speed maximum at the depth of the intrusion noted in Figure 4.4. The 
degree of active mixing which is occurring is illustrated in Figure 4.7. y, the 


temperature variance dissipation rate is defined as 
ee 7 lees, (4.1) 


Where D is the molecular diffusivity of seawater. The maximum in the log of ¥ 
associated with the intrusion at 58 meters indicates an increase in active mixing of 
nearly 4 orders of magnitude near the top of the thermocline as compared to the fossil 
mixed layer. As noted earlier, signatures such as those in Figures 4.4 through 4.7 were 
sporadic, though common, and were often sustained for long periods of time. The 
coincidence of the high x layer and the salinity-compensating temperature inversions 
suggests that diffusive processes are dominating the mixing at these depths. This 
mixing activity will be examined in later studies which emphasize turbulent mixing in 


the upper thermocline. 
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B. HORIZONTAL VARIABILITY 

Figure 4.8 illustrates two of the more salient features of the time series of 
temperature at 10, 20, 40 and 100 meters depth from 6-12 October. Most evident, 1s 
that the upper mixed layer warmed from 6-8 October when incomplete daytime mixing 
of the seasonal mixed layer was observed. During the remainder of the period, the 
upper ocean experienced virtually uninterrupted cooling with well mixed conditions 
prevailing throughout the mixed layer. These heating and cooling trends, which must 
be clearly reflected in an analysis of the total heat content, are examined in the next 
section. Barring significant horizontal advection or vertical temperature redistribution 
into the mixed layer from the thermocline below, the trends must reflect the surface 
forcing outlined in the previous chapter. 

Also apparent in the time series of Figure 4.8 are the horizontal temperature 
trends which exist within the domain and evolve from pattern to pattern. They are 
difficult to quantify by inspection due to the irregular spacing which exists between 
similar measurements at a given depth. That is why the vertical lines which divide the 
time series into individual patterns are not equally spaced. The x-axis is the profile 
number and not distance or time. The variability of the ship speed and the towyo 
winch profile range is responsible for the difficulty in quantifying the trends. The 
character of these gradients or trends is more variable as depth increases below 20 
meters. Inspection of Figure 4.8 clearly shows that there is vertical coherence of the 
time series at 10 and 20 meters depth. As one moves below the active surface nuxing 
layer this coherence decreases. The structure of the time series 1s also clearly different. 
It is tmportant to note that the range of the temperature variability along each leg at 
40 meters is less for the first three patterns of the survey than is observed in the 10 
meter time series for the same patterns. Care must be exercised when comparing the 
time series since each graph 1s scaled between the maximum and minimum temperature 
along the track covered and the ranges differ considerably both in depth and in time. 
Below the base of the mixed layer the presence of internal wave activity prevents the 
simple extraction of horizontal gradients. Correlation of the temperature and salinity 
time series at 10, 20 and 40 meters meters can be seen by comparing Figures 4.8 and 
4.9. 

To understand the temperature trends which Figure 4.8 depicts, time series from 
individual butterfly shaped patterns were plotted with the start and end of each leg of 


the sampling pattern marked. Figure 4.10 depicts similar time series for the first 
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Figure 4.8 10, 20, 40 and 100 meter Temperature, 6-12 October. 
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10, 20, and 40 meter Salinity, 6-12 October. 
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Figure 4.10 10, 40 and 100 mcter Temperature, pattern 4. 
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complete butterfly pattern of the survey. The three vertical lines within the each of the 
graphs divide the patterns into their four component legs. The figure further illustrates 
the degree of spatial variabilitv which exists within the domain and emphasizes that the 
earlier gradients appear stronger at 10 meters than at 40 meters. [t also demonstrates 
in greater detail that the horizontal structure in the upper thermocline is dominated by 
internal wave activity and that mean horizontal gradients which may exist there can't 
be determined in this format. 

Table 6 illustrates the magnitude of the temperature trends and the standard 
deviations of those trends from the starting to the ending profile number along the first 
leg (northwesterly) of each pattern at depths of 10, 40 and 100 meters. The standard 


deviations tabulated with each trend represent the square root of the residual variance 


oie > 
TEMPERATURE TREND (MILLIDEGREES/LEG) , LEG | 


Pattern 10 Meters 40 Meters 100 Meters 
4 -253 2 38 46212 -§ + 145 
5 see a? $3 au -725< 166 
6 -46 + 33 3810 4332 
7 2538 Sec 2043 1) 
8 66217 6814 62324164 
9 iQ e260 pees lo 02 
10 95 = 12 “AQ =e 258 +129 
ie ee ol 33) 399 + 190 
12 S72 A416 95+125 
13 63+9 6110 -302 + 90 
14 Cecelia 40-0 0) lle 
15 ioe 46 2 lee 109 
16 425 a4 2G -985 + 143 


of the leg after the trend is removed. In patterns 4-6 a significant negative temperature 
gradient apparently exists at 10 meters. Since the upper ocean was diluted by rainfall in 
the first several days of the experiment, these gradients are very likely due to the 
advection of warmer, dilute water from the northwest in connection with the prevailing 
north to northwesterly wind induced flow. Gradients, such as those observed, could be 
attributed to transitory synoptic mesoscale fronts, eddies or filaments if a sub-surface 
gradient were also observed. Since the active mixing layer does not extend below 35 
meters (in some cases, and generally less in others) the gradient is probably not the 


result of mesoscale interactions. It is most likely connected with differential surface 
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heating and the advection of dilute surface water from a region of greater precipitation 
to the northwest. The variabilitv in solar heating which was observed due to clouds 
could account for a portion of horizontal gradients but the principal contributor 
appears to be precipitation. The time series of salinity at 10 meters in Figure 4.9 
begins with salinity increasing which corresponds to the northwest to southeast leg of 
pattern 4. 

Dissimilar temperature trends develop at 10 meters in Table 6 after the period of 
precipitation. These gradients are more similar to their counterparts at 40 meters and 
in several cases are comparable. This occured when the convective mixing spanned the 
entire fossil mixed layer from 8-12 October. The higher variability which exists at 100 
meters due to internal wave activity in the upper thermocline is illustrated in column 3 
of Table 6. At 10 and 40 meters internal waves exhibit a minimum of visible influence 
on the thermal structure due to the lack of stratification in the mixed layer. 

The temperature trends at 10 meters along the complete sections of each leg 


containing profiles from near the surface to 120 meters are summarized in Table 7. 
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Complete timeseries exist along legs 1,2 and 4 while only the first half of leg 3 was 
utilizable. The second half of leg three was interrupted with deep tows and tows at 
constant depth and was necessarily omitted. To appreciate the extent of the gradients 
involved, the length of the survey legs described above average 11.2, 8.0, 5.6 and 8.0 
kilometers for legs 1-4 respectively. Patterns 4-6 of Table 7 depict the stronger 
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temperature trends which existed along all four legs during the pre-convective period 
for the mixed layer. In subsequent patterns the large standard deviations indicate that 
significant local temperature gradients often overshadow the smaller trends. 

Table 8 summarizes the time series of temperatures at 40 meters.. The stronger 


trends observed in the first three patterns of Table 7 at 10 meters are smaller in Table 
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8. Although convective conditions extend throughout the entire mixed layer after 
patterns 4-6, the trends at 40 meters do not accurately reflect those at 10 meters for 
several days. Near the strongest convection period toward the end of the study, the 
fossil mixed layer is nearly homogeneous and there is almost a 1:1 correspondence in 
trends at 10 and 40 meters. This can be seen by comparing the trends in Tables 7 and 8 


for patterns 15 and 16 where there 1s very little horizontal structure. 


C. MIXED LAYER HEAT CONTENT 

In the previous chapter on surface fluxes, Figure 3.7 showed two distinct periods. 
From 06 to 08 October, a diurnal cycle occurred in which Q, alternated between +5 
MJ/m? during the day and -5 MJ/m? at night. Over the next four days, 8-12 October, 
there was a net loss of heat at the ocean’s surface of 35 MJ/m?. This cooling can also 
be seen in Figure 4.11, the time series of total heat in the mixed layer. 

Equation 1.2 was integrated to the base of the mixed layer to produce the time 
series in Figure 4.11. To ensure that colder temperatures below the base of the mixed 


layer were not included in the integration, the average temperature of the layer was 
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integrated to an arbitrary depth (1 meter) above the top of the thermocline. This 
method accounts for approximately 2 percent heat content loss over an average depth 
of 50 meters in a completely isothermal mixed layer. The discrepancy is preferable to 
the cold water contamunation that that would occurr across the base of the mixed laver 
due to internal waves if the integration were performed all the way to the mixed layer 
depth. 

Based on a 14 percent maximum error in cumulative net heating, the cumulative 
Qo at the ocean surface was 355 MJ/m*. Integrated over the average depth of the 
mixed layer (50 meters) the surface flux represents a loss of 0.7+0.1 MJ/m?° of total 
heat content. Figure 4.11 shows that during the same period the loss of total heat in 
the mixed laver is approximately 0.6 to 0.8 MJ/m?. This is consistent with the latent 
and sensible heat loss being the primary heat transfer mechanisms during the period. 
The moderate wind conditions during PATCHEX were clearly sufficient enough to 
remove a significant amount of heat from the ocean surface and to mix the warm, 
semi-stratified, daytime surface layer into the entire mixed layer. 

It is important to consider how much of the wind stess translated into vertical 
entrainment of colder water below the base of the mixed layer. Large er. al. (1986) 
determined that strong storms accounted for most of the autuminal cooling of the sea 
surface. During a storm event, one would expect vertical entrainment to be a primary 
factor in the overall heat content of the mixed layer. The period of this study was one 
of moderate winds where low relative humidity coupled with decreased insolation 
allowed substantial cooling of the ocean surface without significant erosion of the 
mixed layer (or entrainment of cooler water). Given the errors in the fluxes, however, it 
is iniportant not to discount the possibility of vertical entrainment at the base of the 
mixed layer. Vertical entrainment might account for the observed error in the balance 
or it may be offset by horizontal advection. Barring the presence of mesoscale currents 
to provide net transport, horizontal advection is a questionable cooling mechanism for 
this period. In summary, since there were no storm episodes during PATCHEX, it 1s 


understandable that the heat content balances well with the large surface fluxes. 
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Figure 4.11 Mixed layer heat content, 6-12 October. 
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V. HORIZONTAL GRADIENTS BY OBJECTIVE ANALYSIS 


A. THEORY AND TESTING OF THE O-A 

The space-time series of temperature at constant depth were mapped on to a 
uniform grid using an objective analysis routine written by J. Smith at the Naval 
Postgraduate School and recently adapted for PATCHEX. This analysis was 
undertaken in order to develop quasi-synoptic snapshots for each of the 13 consecutive 
Survey patterns this thesis employs. The algorithm 1s based on the work of Bretherton 
et al (1976) which used MODE (Mid Ocean Dynamics Experiment) data to develop a 
successful oceanographic application of interpolator/extrapolator techniques which had 
been developed earlier for meteorological analyses (Gandin, 1965). The program 
employs a simple 3 dimensional radial search in X and Y space and in time to find 
seven scalar values closest to each gridpoint. Values are taken from each profile at a 
given depth, y,, equal to the magnitude of the scalar of interest (temperature in this 
case) and added to an estimate of the random noise associated with the measurement, 


fe, 
ee On c.. (&5.1) 


The objective analysis then uses those values to estimate the value of the scalar and its 
error at successive grid points within the PATCHEX domain. To produce this spatial 
depiction, the algorithm uses pre-determined correlation scales to map the scalar onto 
a uniform grid. 

Determining the correlation length for the data in each pattern was complicated 
by the irregular spacing and sampling interval between successive measurements of 
scalar quantities at a given depth. To determine the autocorrelation length of the data, 
an autocorelation function was generated by calculating a radial lag between successive 
data points along the survey track and every other point in the data field. The 


autocorrelation was then calculated using the standard formula 
Ry = (%- Wp br > WS" (5.2) 
where r 1s the radial lag. This method assumes an isotropic data field. 
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A graphic depiction of the radial correlation function was produced by dividing a 
8 kilometer radius into 20 bins and accumulating autocorrelation values into respective 
bins every 0.4 kilometers from the origin. Figure 5.1 depicts the correlation function 


for a typical survey pattern. The modal value of the correlation radius for the 13 
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Figure 5.1 Typical correlation function, pattern #4. 


patterns that this study employs is approximately 6 kilometers. Since the data set is 
finite in size and the algorithm divides the abcissa of 8 kilometers into 20 bins, the 
method docs not perfectly reproduce the shape of a continuous autocorrelation curve. 
It does, however, effectively determine the radius at which the data becomes 
decorrelated. 

A default correlation length, Ry, of 12 kilometers was applicd in the objective 
analysis. The correlation function it uses is double sided and our correlation radius as 


determined earlier is 6 kilometers. Specifically, the correlation function used in the 
objective analysis is 


F = (1-R) eRe) (5.3) 


Where the three dimensional correlation scale, Ry is defined as 
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Ree aoe) aye 17/72, (5.4) 


Note that the correlation function includes a time scale. In this analvsis, it was 
defeated, assuming that each pattern was sampled synoptically. The average elapsed 
time of each survey pattern ranged between 6-11 hours. § hours was roughly the 
average. Earlier patterns were surveyed more slowly than later ones. 

After the data were run through a convolution filter, the default correlation scale 
of 12 kilometers was applied along with an error coefficient, €. = 0.2. This value of €. 
represents the average percentage of noise along each leg of each pattern in the survey. 
Testing showed that the objective analysis is not particularly sensitive to €. < 0.5. The 
average of the high end of the observed noise (20 percent) proved adequate in 
reproducing the resolvable temperature trends in each pattern without over-smoothing. 

Before the actual data could be mapped, however, additional testing was 
undertaken. The O-A routine was tested with 2 fields, a south to north temperature 
gradient and a temperature spike imbedded in the gradient. The spike proved useful in 
tuning the filtering and correlation coefficients for later application in the actual 
analysis. It was located in the southwest corner of an 18.0-18.4°C south to north 
temperature gradient. The V,T of 0.4 degrees over 12 kilometers was a representative 
temperature gradient in the same range as the time series observed in PATCHEX. The 
spike consisted of four adjacent points having a value of 18.3°C. The overall size of 
the input grid was 12 by 12 kilometers and the data were 400 points spaced 0.6 
kilometers apart. The input grid is contoured in Figure 5.2. It depicts a maximum 
closed temperature contour of 18.24°C surrounding the spike. The S-N temperature 
gradient is also apparent in the figure. 

The first test of the objective analysis was made using a correlation length, Ro, of 
2 kilometers to more easily resolve the 0.6 kilometer spaced data. Results of that O-A 
Were then contoured on the same grid used for the input field and are depicted in 
Figure 5.3. The input array had been subjected to various filter and error coefficents as 
well as different correlation lengths. The output array was best resolved using 
atriangular convolution filter over £4 points, an error coefficient of 20 percent and a 2 
kilometer correlation length. Using these input parameters, an 18.20 maximum closed 
contour was analyzed.This represents an an 8 percent error in the estimate of the spike 
test function by objective analysis. Although various combinations of the input 


parameters were tested to produce the final field, the figure represents the best mix of 
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Figure 5.2 12kmN X 12 km E, O-A test input field. 


those parameters which are still germane to PATCIIEX. The convolution filter was 
used to reduce high frequency internal wave noise in the input field. The error 
coeflicient was estimated from the high frequency content of each PATCITEX pattern 
used. The 2 kilometer correlation length tested the ability of the algorithm to better 
resolve localized features. Since the actual correlation length scale of the PATCHIEX 


data is 12 kilometers, the data were re-tested using the higher value. This resulted in an 
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Figure 5.3 12 kmN X 12 kmE, O-A test output field, ae 2 km. 


expected loss of resolution of higher spatial frequency features. Since the actual scale of 
the test field was less than a kilometer, a 12 kilometer correlation scale is inappropriate 
for the high spatial frequency of the test data. The results of the analysis can be seen in 
Figure 5.4 where the largest resolved closed contour ts 20.16°C. The artifact at the 
center of the contour indicates further error in the analysis. The next section will show 


that a larger R, is more applicable to the PATCIIEX field, which is well correlated to 


12 kilometers. 
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Figure 5.4 12 kmN X 12 kmE, O-A test output field, Ry = 12 km. 
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B. RESULTS OF THE OBJECTIVE ANALYSIS OF PATCHEX DATA 


The objective analvsis was applied to the PATCHEX data using the input 
parameters described in the last section. The correlation scale was 12 kilometers. The 


results of the analysis are depicted in Figure 5.5. Also overlaved in the figure is the 





< 
~ 
Ww 
© 
= 
Parse vente eCRaGURPE ORG Ce). Box 4 12M} Relss he Kin 
ngs lovese; (lax: 20515 Contaur interval: ,.G2Prrg 
Fast ————— 
Figure 5.5 


10 meter temperature, pattern 4, 14 km N X 14 kmE, Ry = 12 km. 


butterfly shaped ship’s track. The objective analysis was mapped on to a 14 X 14 
kilometer grid of 1600 (40 X 40) points. The error field was also computed and 1s 


63 


depicted in Figure 5.6. It indicates that the objective analysis was able to produce 
accurate estimates of the temperature along the survey track but was only marginallv 
successful in calculating the rest of the field. As would be expected, the regions outside 
the track and at the center of the two halves of the pattern exhibited especially high 


errors. 
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Figure 5.6 10 meter temperature error, pattern 4, 14 kmN X 14 kmE, Ry = 12 km. 
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Inspection of the analysis in Figure 5.5 shows that the along-track temperature 
trends correlate extremely well with those in Table 6 of Chapter IV. It further shows 
that the general horizontal temperature gradient is decreasing from the northwest to 
the southeast which also appeared in Table 6. Warmer, more dilute surface water 
clearly exists in the northwestern regions of the domain as was suspected in the earlier 
description of the raw temperature time series. The most erroneous feature in the 
objective analysis is the extreme gradient located at the center of the eastern half of the 
survey pattern. This is substantiated in the error analysis of Figure 5.6. 

Various attempts were made to adjust the input coefficients to better resolve the 
fields by objective analysis. Attempts were even made to decrease the correlation time 
scales and reduce the effect of temporally distant data. This approach was rejected after 
the incomplete fields provided extremely erroneous and noisv results. The effect of 
decreasing the correlation length in order to better resolve localized trends are seen in 
Figures 5.7 and 5.8. Using a 6 kilometer correlation length, Figure 5.7 was able to 
better isolate the trends along the survey track but increased the off-track noise. Also, 
the actual values of the along-track trends resulting from applying a decreased 
correlation length were no more accurate than those in Figure 5.5. Further reduction of 
the correlation length to 3 kilometers only exacerbated the error away from the track. 
This is depicted in Figure 5.8. 

Having confined our choice of correlation length, filter coefficient and error 
criteria, it was decided to apply the OA only to depict the trends along those survey 
patterns whose trends were substantially greater than their associated noise values 
(variances) as tabulated in Tables 7 and 8. This decision proved correct as the OA had 
performed poorly in analyzing pattern 4 which contained the most significant 
temperature trends. Testing of subsequent patterns indicated that as the mixed layer 
became more homogeneous, with less defined large-scale gradients, the OA lost the 
ability to discriminate the trends. Since the first three patterns of the survey covered a 
period when a surface layer with significant horizontal gradients existed, those patterns 
contained trends greater than the order of the high frequency noise within the domain. 
Pattern 6 contained the last clearly resolvable horizontal trend at 10 meters. 
Following pattern 6, the trends are lost in the noise and the O-A begins to create 
artifacts. Despite the noisier subsequent patterns, Figures 5.9 and 5.10 depict the 
decreasing temperature trend from the northwest to the southwest corner as in Figure 
a 
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Figure 5.7 10 meter temperature, pattern 4, 14 km N X 14 kmE, Rp = 6 km. 
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Figure 5.8 10 meter temperature, pattern 4, 14 km N X l4 km E, Ry 
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10 meter temperature, pattern 5, 14 km N X 14 kmE, Ry = 12 km. 
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Figure 5.10 10 meter temperature, pattern 6,14 kmN X l4 km E, Ry = 12 km. 
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Similar results were obtained using the data at 40 meters. Again, the first three 
patterns of the study provided the best results. As in the 10 meter analyses, the 
gradients along the track were mapped correctly but the rest of the domain suffered 
very high errors. The results of the analyses at 40 meters for the first three patterns 
can be seen in Figures 5.11 through 5.13. Note that these analyses resolved the overall 
trends at 40 meters more clearly in Figures 5.12 and 5.13 than in Figure 5.11. Since this 
region is below the actively mixing surface layer and near the base of the fossil mixed 
layer, communication with the upper thermocline may temporally alter the temperature 
field. Internal waves, convective compensation processes near an intrusion and vertical 
transport over the 8 hour survey cycle can all contribute to quickly change the 
horizontal temperature distribution near the base of the mixed layer. Despite these 
variations, which are extremely small in magnitude and not easily resolved, the O-A 
was still able to portray the general horizontal distribution of temperature at 40 meters 
in the first 3 patterns. 

Inspection indicates that there is a positive west to east temperature gradient in 
each of these patterns. This further implies that the surface layer and the fossil mixed 
layer are completely decoupled during the first two days of the study. The temperatures 
at 40 meters are also cooler than the temperatures at 10 meters during the two day 
period when convective mixing of the fossil mixed layer and the surface layer was 
curtailed. These results support the hypothesis that a warmer, saltier, northwesterly 
moving layer of surface water inhabits the PATCHEX domain on 6-8 October 1986. 
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40 meter temperature, pattern 4, 14 km N X 14 kmE, Ry = 12 km. 
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Figure 5.12 40 meter temperature, pattern 5, 14 km NX l4 kmE, Ry = 12 km. 
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Figure 5.13 40 meter temperature, pattern 6, 14 kmN X 14 km E, Ry = [2 km. 
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VI. SUMMARY AND RECOMMENDATIONS 


The period from 6-12 October 1986, was examined as a subset of the PATCHEX 
experiment. During this time, there were two distinctly different periods of heat flux 
and horizontal temperature structure in the mixed layer. To understand the differences 
between the 6-8 October time-frame and the four days that followed, the surface fluxes, 
heat content, and horizontally-varving, constant-depth, temperature time series were 
examuned. 

The arrival of the R/V POINT SUR was at the end of a period of precipitation. 
The influence of an intensifying anticyclone sets the stage for the conditions that 
followed. The wind stress averaged 0.1 N/m? until the last day of the study. During the 
first two days, a -150 to 500 W/m? diurnal fluctuation existed in the net surface heat 
flux. An increasing dewpoint depression over the next 4 days was responsible for 
periods of over 100 W/m? of latent heat loss. The maximum daytime insolation 
dropped from approximately 600 MJ/m* on 6 October abruptly to 150 MJ/m? on 9 
and 10 October. This kept the surface from regaining the large amount of energy it had 
lost to latent and sensible heating. Cumulative time series showed a diurnal fluction of 
+5 MJ/m? early in the study and a net loss of 35 + 5 MJ/m? by the end of the study. 

A detailed vertical and horizontal thermohaline analysis of the mixed layer was 
undertaken. It suggested that from 06 to 08 October, warmer (0.1°C to ).2°C) and less 
saline (0.04 to 0.16 parts per thousand) water comprised a 10-35 meter surface slab 
over a fossil mixed layer extending to 42-62 meters depth. During the second part of 
the study, stirring from a 0.1 N/m* mean wind stress, and convective muxing from the 
large surface heat flux maintained the vertical homogeneity of the mixed laver. The 
analysis also indicated that the upper thermocline was active. Strong stratification, 
internal waves, thermohaline intrusions, active mixing, and a 0.3 knot current, several 
meters below the base of the mixed laver, were noted throughout the study. 

One of the most important results of this study was the discovery of horizontal 
temperature gradients in the time series of temperature at 10 and 40 meters. These 
gradients at 10 and 40 meters, were overlaid upon a 2-day warming trend followed by 4 
davs of cooling. These longer term heating and cooling trends, as opposed to the short 


term trends in the time series, were reflected in the surface fluxes. The horizontal 
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gradients were strongest from 6-8 October and were different at 10 and 40 meters. At 
1Q meters they were strongest and increased most rapidly from the northwest to the 
southeast with decreasing temperature. At 40 meters the gradients were weaker. Also, 
the increasing of the gradient to the southeast was not observed. After 8 October, 
convective overturning of the layer reduced the magnitude of the gradients to the order 
of the temporal noise. 

The time series of heat content in the mixed layer showed the same trends as the 
time series of surface heating and fixed depth temperature over the two periods of the 
study. During the first period, the heat content increased by approximately 0.5 MJ/m?, 
and peaked ober, on 8 October. It then decreased by 0.6-0.8 MJ/m? by 12 October. 
This large heat content loss in the mixed layer balanced well against a 0.7 + 0.1 
MJ/m> surface loss, however, the importance of horizontal advection and vertical 
redistribution from entrainment, were not discounted. 

Objective analysis of the 10 meter temperature timeseries on 6 October indicated 
that a temperature gradient in excess of -0.25 m°C/km existed along the northwest to 
southeast transect of the area. It persisted, but weakened, in two analyses of later 
transects. At 40 meters, the gradients were weaker and positive from west to east. After 
the first two davs of the study, the mixed layer was virtually isothermal, and the weak 
gradients were indiscernable from the high frequency noise in the time series at 10 and 
40 meters. 

After examining the persistence of temperature gradients in the time series and 
their further resolution by objective analysis and coupling that information with the 
vertical thermohaline character of the mixed layer, it appears that from 6-8 October 
there was a two layer mixed layer. Wind stress and extreme heat loss mixed those 
layers to a virtually isothermal and isohaline state beyond 8 October. The data suggest 
a rapidly adjusting density interface originally caused by precipitation and the heating 
of surface water which 1s eliminated by wind stress and convective mixing. The study 
neither confirms the existence of a mesoscale feature nor does it dtscount their 
existence. The slope of the maximum temperature gradient (25 m°C/km) exceeds that 
Which would be observed in conjunction with a mesoscale feature. It 1s more likely the 
result of density compensation at the interface between a pool of fresh surface water 
and the fossil mixed layer. The existence of that extreme temperature gradient, 
however, can easily mask the presence of a mesoscale signature. To clarify this matter, 


It 1s proposed that the Expendable Current Profiler data, large scale CTD data and 


fe 


shipboard ADVP be examined to discern the large scale interactions surrounding 
PATCHES 

Once the mean velocity wfeldvis obtained frome FLIP, and whenmeic wertical 
entrainment velocity can be estimated for the R/V POINT SUR data, a more complete 
closure of the heat balance in the upper ocean should be attempted. Also, to determine 
the nature of the internal wave field from the PATCHEX data set (both in space and 


time) a maximum entropy spectral analysis might be undertaken. 
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